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Abstract 

Ice-stream  discharge  fluctuations  constitute  an  independent  means 
of  forcing  unsteady  ice-shelf  behavior,  and  their  effect  must  be 
distinguished  from  those  of  oceanic  and  atmospheric  climate  to 
understand  ice-shelf  change.  Ice-stream-generated  fluctuations  of  a 
one-dimensional  ice  shelf  are  found  hereto  propagate  through  the  ice- 
shelf  environment  along  the  two  characteristics  of  the  hyperbolic 
governing  equations.  One  characteristic  permits  instantaneous 
transmission  of  grounding-line  velocity  changes  to  all  points 
downstream.  The  other  characteristic  represents  slow  transmission  of 
grounding-line  thickness  changes  along  Lagrangian  particle  paths. 
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Introduction 

Unsteady  thickness  and  flow  observed  on  the  Ross  Ice  Shelf  are 
thought  to  result  primarily  from  fluctuations  of  ice  streams  and  outlet 
glaciers  that  feed  the  ice  shelf  at  its  inland  boundaries. (Shabtaie  and 
Bentley,  1987;  MacAyeal  and  others,  in  press).  The  initial  effects  of 
climatic  change  thus  must  be  discriminated  from  an  ever-present 
background  of  ice-stream-generated  variability.  To  characterize  this 
variability,  we  determine  the  general  characteristics  of  ice-thickness 
and  flow  anomalies  produced  by  time-varying  grounding-line  discharge.  We 
restrict  our  examination  to  anomalies  of  an  ideal,  floating  ice  tongue 
to  facilitate  and  analytic  treatment  of  the  governing  equations.  Our 
results  provide:  (i)  a  conceptual  basis  for  modeling  thickness  and  flow 
variations  of  West  Antarctic  ice  shelves  (such  as  developed  in  a 
companion  paper,  MacAyeal  and  Lange,  unpublished) ,  and  (ii)  an 
understanding  of  how  ini  and- ice  discharge  may  compete  with  oceanic  and 
atmospheric  conditions  to  force  such  variations. 


Characteristic  Equations 

Thickness  H(x,t)  and  longitudinal  velocity  U(x,t)  for  the  ideal 
ice-shelf  geometry  (x>0)  shown  in  Figure  1  satisfy 

* 

atH  +  UaxH  -  M  -  A"Hn+l  (1) 

axU  -  A"Hn  (2) 

subject  to  boundary  conditions  (specifying  H  and  u  as  functions  of  time) 

>ri  F ol* 

at  the  grounding  line  x-0  (note,  3  denotes  the  x-derivative,  etc.). 

*  A&I 

Equations  1  and  2  are  derived  in  the  Appendix  and  stem  from  mass- 

v' .  a 

continuity  and  stress-equilibrium  principles,  subject  to  certain 
simplifying  assumptions.  Dimensionless  variables  are  defined  in  Table  I 
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using  scales  typical  of  Antarctic  conditions  (Holdsworth,  1985).  iiuy 
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Steady-state  solutions  to  Equations  1  and  2,  Hs(x)  and  U  (x), 
satisfying  Hs(0)=l  and  Us(0)=l  at  the  grounding  line  were  found  by  Van 
der  Veen  (1983)  and  are  written  (for  M>0): 

Hs(x)  -  [(l-7)/(l+Mx)n+1  +  (3) 

Us(x)  =  [1  -  7  +  7(l+Mx)n+1](1/(n+l!)  (4) 

where  7=An/M.  Hs(x)  monotonically  decreases  with  increasing  distance 
from  the  grounding  line  and  approaches  a  constant  value  7'1/(n+1)  as  x-®. 
This  asymptotic  thickness  (approximately  200  m  for  typical  Antarctic 
conditions)  is  achieved  because  ice-shelf  thinning  by  horizontal 
spreading  is  balanced  by  snow  accumulation  M.  Us(x)  monotonically 
increases  with  distance  from  the  grounding  line,  but  does  not  approach 
an  asymptotic  limit  as  x-*«. 

Deviations  from  steady  state  caused  by  fluctuations  in  H  and  u  at 
the  grounding  line  can  be  computed  directly  from  Equations  1  and  2. 
Alternative  forms  of  Equations  1  and  2,  called  the  characteristic 
equations,  are  useful,  however,  for  displaying  fundamental  properties  of 
the  solution.  They  are  derived  by  defining  new  spatial  and  temporal 
coordinates  that  produce  ordinary  differential  equations  for  thickness 
and  velocity  (the  method  of  characteristics)  and  are  written: 


aTx  -  udT t  -  o 

(5) 

d(t  -  0 

(6) 

aTw  -  (M-An/r+1)aTt  -  o 

(7) 

-  a n«"afx  -  o 

(8) 

where  W(£,r)  and  (/(£,r)  represent  thickness  and  velocity,  respectively, 
as  functions  of  alternative  coordinates  £  and  r  that  satisfy 


age 

The  curves  £  =  £  and  r  =  ro  in  the  x,t  plane  represent  characteristics 
of  the  ice-shelf  system  along  which  changing  conditions  propagate.  These 
curves  can  be  determined  by  solving  Equations  5  and  6  for  x(£,t)  and 
t(r)  (t  and  r  can  be  taken  as  identical  in  the  present  application;  this 

gives  drt  =1).  As  shown  by  Equation  5,  the  characteristic  £  =  £ 
represents  the  trajectory  of  a  particular  ice  column  (which  may  be 

labeled  by  its  ^-coordinate  £o)  in  the  x,t  plane.  The  thickness  of  this 

particular  ice  column,  tf(£=£Q,T),  is  determined  as  a  function  of  t  by 

Equation  7.  The  velocity  of  this  ice  column,  (/(£=£0,r),  is  found  by 

integrating  Equation  8  over  the  domain  where  ^(r)  is  the  £- 

label  of  the  ice  column  currently  at  the  grounding  line.  This  interval, 

in  x-coordinates,  spans  the  distance  between  the  grounding  line  and  x  = 

x(£0,r). 

The  characteristic  equations  illuminate  several  basic  aspects  of 
time-dependent  ice-shelf  adjustment.  Equation  7  indicates  that  the 
thickness  of  a  given  ice  column  (^-constant)  evolves  independently  of 
the  ice  thickness  elsewhere.  This  independence  is  an  aspect  of 
unconfined  ice  shelves  only.  In  circumstances  where  coastal  confinement 
introduces  resistance  to  seaward  flow,  evolution  of  a  given  ice  column 
depends  on  thickness  conditions  downstream. 

Another  basic  aspect  of  ice-shelf  behavior  illuminated  by 
Equations  7  and  8  is  that  deviations  from  steady-state  are  limited  in 
space  and  time  by  the  appropriate  characteristic  curves.  If,  for 
example,  H  or  U  at  x»0  are  impulsively  changed  at  t*0  and  the  initial 
conditions  are  in  steady  state,  unsteady  conditions  will  prevail  only 
within  the  region  of  the  x,t  plane  bounded  ty  the  curves  ^-constant  and 
T*constant  that  emanate  from  x-t-0.  This  partition  of  the  x,t  plane 
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into  steady  and  unsteady  regions  is  shown  schematically  in  Figure  2. 
Change  propagates  downstream  instantaneously  along  the  r=constant  curve 
because  velocity  at  any  point  is  determined  by  spatial  integration  of 
the  instantaneous  ice-thickness  distribution.  Equilibration  with  new 
forcing  conditions  is  completed,  however,  in  a  delayed  fashion  following 
the  ^constant  curve.  Steady-state  conditions  prevail  in  the  wake  of  the 
ice  column  that  was  at  the  grounding  line  at  the  instant  t=0.  This  ice 
column  follows  the  path  determined  by  the  ^constant  characteristic. 
Transient  Ice-Shelf  Profiles 

We  demonstrate  basic  ice-shelf  response  by  computing  ice  thickness 
and  velocity  anomalies  generated  by  prescribed  scenarios  of  unsteady 
grounding-line  discharge.  These  anomalies  are  defined  by 

h (x , t )  -  H(x,t)  -  Hs(x)  (11) 

u(x,t)  -  U(x,t)  -  US(X)  (12) 

where  Hs(x)  and  Us(x)  are  the  steady-state  thickness  and  velocity 
profiles  defined  by  Equations  3  and  4,  respectively.  To  simplify  the 
calculations,  we  require  h  and  u  to  be  small  and  use  a  linear  expansion 
of  Equations  1  and  2: 

3th  +  Us3xh  -  -(n+l)AnHsnh  -  3xHsu  (13) 

3xu  -  nAnHs"-lh  (14) 

With  this  simplification,  the  two  generic  characteristic  curves 
E-constant  and  r-constant  associated  with  steady-state  conditions  may  be 
applied  as  approximations  to  the  characteristic  curves  associated  with 
unsteady  conditions. 

Boundary  conditions  at  the  grounding  line  (x*0)  and  initial 
conditions  are  specified  to  represent  three  basic  scenarios  of  ice- 
stream  discharge  fluctuations  (amplitudes  and  time  constants  are 


arbitrary  for  demonstration): 
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scenario  I  -  u(x,0)  =  h(x,0)  -  0.0  (15) 

u(0,t)  =  h(0,t)  =  0.1  0<t<0.2 

0.0  t>0.2 

scenario  II  -  h(0,t)  =  0.1  sin(ut)  (16) 

u(0,t)  -  0.0 

scenario  III  -  h(0,t)  -  0.0  (17) 

u(0,t)  =  0.1  sin(wt) 

where  w«2*/0.5  .  Scenario  I  represents  an  episode  of  increased 
grounding-line  thickness  and  velocity  that  is  sustained  for  a  brief  time 
interval.  Scenarios  II  and  III  represent  periodic  fluctuations  in  h  and 
u  (we  do  not  consider  a  variety  of  forcing  frequencies  in  the  present 
study;  see,  however,  MacAyeal  and  Lange,  unpublished).  Initial 
conditions  for  scenarios  II  and  III  are  not  required  because  the 
grounding  line  forcing  is  assumed  to  proceed  ad  infinitum  (in  practice, 
we  apply  an  initial  condition  of  h(x,0)-0  and  u(x,0)»0,  and  proceed  with 
time  integration  until  all  manifestations  of  the  initial  conditions  have 
dissipated) . 

Finite-difference  solutions  to  Equations  13  and  14  under  the  three 
forcing  scenarios  described  above  were  produced  on  the  interval  0<x^l 
and  over  a  sufficient  time  interval  to  display  salient  features  (and  to 
dissipate  transients  induced  by  initial  conditions  in  periodic  forcing 
scenarios).  The  finite-difference  form  of  Equation  13  possessed 
implicit  time  steps  (At-0.01)  and  centered  space  derivatives  (Ax-0.01), 
except  at  x-1.0  where  an  upwind  space  derivative  was  used.  Equation  14 
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was  integrated  at  each  time  step  over  x  using  centered  space  derivatives 
(except  upwind  at  x-l.O). 

Salient  features  of  scenario  I  are  displayed  by  the  contour  maps 
of  thickness  anomaly  h  and  volume  flux  anomaly  q,  defined  by 

q(x,t)  -  u(x,t)Hs(x)  +  h(x,y,t)Us(x)  (18) 

in  Figure  3.  In  Figure  4,  the  x,t  plane  on  which  these  maps  are 
presented  is  divided  into  five  regions  by  the  four  characteristic  curves 
that  emanate  from  the  grounding  line  at  the  start  and  end  of  the 
modified  discharge  episode.  In  regions  I,  III  and  V  shown  in  Figure  3, 
ice-shelf  conditions  are  in  steady  state  with  current  discharge 
thickness  and  velocity  [Hs(0)  +  h(0,t)  and  Us(0)  +u(0,t)].  In  regions 

II  and  IV,  conditions  are  unsteady.  Comparison  of  the  maps  of  h  and  q 
indicates  that  thickness  anomalies  are  organized  primarily  along 
^-constant  trajectories  and  that  flux  anomalies  tend  to  display  rapid, 
discontinuous  changes  along  r*constant  trajectories. 

Ice-thickness  anomalies  generated  in  response  to  scenarios  II  and 

III  (periodic  forcing)  indicate  two  contrasting  patterns  shown  in  Figure 
5.  For  oscillatory  grounding-line  thickness  (scenario  II),  h  is 
maximum  at  the  grounding  line  and  decay  monotonically  downstream. 

Ridges  and  troughs  of  the  h(x,t)  contours  in  the  x,t  plane  for  this 
situation  extend  along  ^-constant  trajectories.  In  contrast,  the 
maximum  of  h  generated  by  scenario  III  (oscillatory  grounding  line 
velocity)  occurs  downstream  of  the  grounding  line  at  a  location  that 
depends  on  the  frequency  of  forcing  (the  lower  the  frequency,  the 
farther  downstream) . 

Comparison  of  scenarios  II  and  III  suggests  that  phase 


relationships  between  thickness  and  velocity  oscillations  of  a  given  ice 
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stream  Is  critical  in  determining  the  location  of  maximum  ice  thickness 
anomalies  on  the  ice  shelf.  This  aspect  of  ice-stream  dynamics  may  thus 
determine  where  ice  rumples  and  rises  are  likely  to  form  in  response  to 
ice-stream  changes. 

Conclusion 

'i 

Discrimination  between  ice-stream  and  climatically  forced  ice- 
shelf  transience  can  be  accomplished  through  the  relationship  between 
ice-shelf  anomalies  and  the  characteristic  trajectories  that  determine 
propagation  of  ice-shelf  change.  Variability  forced  by  the  atmosphere  or 
ocean  would  not  necessarily  be  constrained  to  propagate  along  such 
characteristics,  but  may  propagate  through  the  ice-shelf  environment 
along  pathways  determined  by  oceanic  and  atmospheric  dynamics. 

The  two  types  of  characteristics  described  represent  fast  and  slow 
paths  by  which  grounding-line  influences  are  transmitted  to  the  ice 
shelf  downstream.  Velocity  and  mass-flux  changes  are  propagated 
instantaneously  via  the  fast  path  because  the  velocity  profile  is 
determined  by  an  elliptic  boundary-value  operator.  Thickness  changes 
are  propagated  advectively  via  the  slow  path  because  ice-column 
thickness  evolves  (in  an  unconfined  ice  shelf)  independently  of  ice- 
column  conditions  elsewhere. 

This  research  was  supported  by  the  National  Science  Foundation 
(DPP  85  09451)  and  the  Office  of  Naval  Research  under  Grant  N00014-86-K- 
0034.  We  thank  L.  Morland  for  productive  discussions  at  the  University 
of  Chicago  in  December,  1986. 
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Appendix:  Derivation  of  Governing  Equations 

Mass  balance,  stress  equilibrium,  heat-flow  continuity  and  the 
material  properties  of  ice  govern  the  evolution  of  large  ice  masses. 

For  ice  shelves  that  float  on  an  effectively  frictionless  ocean,  and 
that  span  a  horizontal  distance  greater  than  their  thickness,  these 
principles  are  considerably  simplified  in  that  vertical  variation  of  the 
horizontal  flow  is  negligible.  This  simplification  has  been 
demonstrated  by  a  number  of  means  (Sanderson  and  Doake,  1979;  Mori  and, 
1987;  and  Muszynski  and  Birchfield,  1987),  and  is  generally  supported  by 
observation.  In  the  present  study,  we  justify  this  simplification  and 
derive  the  governing  equations  for  the  ideal  ice-shelf  configuration  we 
consider  by  expanding  the  stress-equilibrium  equations  and  associated 
boundary  conditions  using  52  (the  square  of  the  aspect  ratio).  We 
conduct  this  derivation  under  the  assumption  that  temperature-dependent 
rheological  parameters  may  be  replaced  with  similar  depth-dependent 
parameters  to  approximate  the  prevalent  vertical  temperature  gradient  in 
Antarctic  ice  shelves.  More  precise  treatment  of  thermo-mechanical 
coupling  (such  as  demonstrated  by  Williams  and  Hutter,  1983;  Shumskiy 
and  Krass,  1979;  or  Morland  and  Shoemaker,  1982)  is  not  necessary  to 
achieve  the  limited  objectives  of  this  study.  We  additionally  restrict 
our  analysis  to  two  dimensions  (vertical  and  horizontal  along  the 
longitudinal  axis). 

Stress  equilibrium  in  the  horizontal  (x)  and  vertical  (z) 
directions,  and  dynamic  boundary  conditions  at  the  ice  shelf  surface 
(zs(x),  stress  free)  and  base  (zb(x),  sea-water  pressure  only),  are 
written: 
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-raxp  +  ax(vexx)  +  r2az(  i/ej  =  o 

Zb<Z<Zs 

(al) 

-ra2p  +  ax(vexz)  +  az(vezz)  -  r  -  o 

Zb<Z<Zs 

(a2) 

-axzs(1/exx  -  rp)  +  =  0 

Z=Zs 

(a3) 

-axZs1/exz  +  UQzz  -  rp  -  o 

Z=Zs 

(a4) 

-  rp)  -  5-2vexz  -  (pjps) rax(ub2) 

Z*Zb 

(a5) 

^xZbVQxz  -  "ezz  +  rP  "  -(^i)rzb 

Z*Zb 

(a6) 

where  we  have  replaced  reference  to  the  deviatoric  stress  using 

V  *  2vejj  (a7) 

and  the  effective  viscosity  u  is  defined  so  as  to  represent  Glen's  flow 
law 

v  -  B(z)/[(%(exx2  +  e„2  +  (a8) 

All  variables  are  defined  in  Table  I  and  are  expressed  in  non- 
dimensional  form  (Hutter,  1983).  As  stated  above,  the  non-dimensional 
flow-law  parameter  B(z)  is  assumed  z-dependent,  rather  than  temperature 
dependent,  to  avoid  the  need  to  treat  thermo-mechanical  coupling. 

Systematic  approximations  to  Equations  al  -  a6  are  developed  by 
expanding  all  variables  in  power  series  using  the  small  parameter  &z. 
The  zero-order  approximations  to  Equations  al  -  a6  are 

5z(,/(0)exz{0))  ■  0  7 


zb""<z<zs""  (a9) 


page  11 


y(0)g  (o)  s  Q 


z=zs(0) 


(all) 


-3  7  t°)«t°)e  (0)  +  u(0)e  (0)  .  rp(0)  .  q 

x  s  xz  zz 


z=zs'°> 


(al2) 


j;(0)g  (0)  =  0 


z=z, 


(0) 


(al3) 


-  t'to,e21(0)  ♦  rp(t»  -  -(.fiJeJTz™  Z-ZB(<»  (al4) 


The  zero-order  solution  can  be  obtained  immediately  by  integrating 
Equations  a9  and  alO  subject  to  the  boundary  conditions  all  -  al4: 


e  (0)  =  0 

xz 


p(0)  -  (zs<°>  -  Z)  +  (y(0)/r)  eJ0) 


zb(0)  »  pAp,  -  pJ  z 


(0) 


(al5) 

(al6) 

(al7) 


The  above  solutions  indicate  that:  (i)  vertical  shear  in  the  ice  column 
is  negligible,  (ii)  pressure  increases  linearly  with  depth,  and  (iii) 
the  ice  shelf  floats  in  local  hydrostatic  equilibrium. 

A  subset  of  the  first-order  approximations  to  Equations  al  -  a6 
(for  the  x-component  of  forces)  describes  how  the  zero-order  strain  rate 
eXx(0)=-e2z<°>  varies  with  H(0)-zs(0)-zbt0): 


a  „(o)e  (i)  =  ra  p{0)  -  diA]eJ0)) 


zb(0)<x<zs(0)  (al8) 


axzs(0)(^(0)exx(0)  -  rp(0)) 


(al9) 


w(Q)e  U) 
exz 

w(0)p  U) 

exz 


Z=Z 


(0, 


Mb 


(0)(i/(0)exx(0)  -  rp(0)  -  pyp,rzb(0))  z=zt 


(0) 


(a20) 


Integration  of  Equation  al8  over  z  using  boundary  conditions  a!9  and  a20 
followed  by  integration  of  the  result  over  x  (assuming  that  exx(0)-0  as 
H(0)-0,  and  that  T,pw  and  ps  are  independent  of  x)  gives 


pz  (0)^(0)  =  ^  r  (i-p./pj  H«°J 

where 


v*  to)  = 


B(z')  dz'  /  (ej0)  1‘1/n) 


Substitution  of  Equation  a22  into  a21  gives 


(a21) 


(a22) 


exx(0)  3  M  3  A"H"  (a24) 

where 

A  -  is  r/B*  (1  -  p./pj  (a25) 


is  a  non-dimensional  parameter  that  measures  the  ratio  of  the 
gravitationally  driven  spreading  stress  to  the  stress  required  to  deform 
the  ice  shelf  at  the  chosen  strain-rate  scale  Ug/Lg. 
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Table  1  -  Variables  and  scales 


efinltion 


Ice  thickness 
horizontal  velocity 
thickness  anomaly 
velocity  anomaly 
characteristic  coordinate 
characteristic  coordinate 
horizontal  coordinate 
vertical  coordinate 
aspect  ratio 
time 

accumulation  rate 
ice  density 
gravity 

sea-water  density 
flow-law  constant 
flow- law  exponent 


stress  ratio 


pressure 

longitudinal  strain  rate 
vertical  strain  rate 
vertical  shear  rate 
effective  viscosity 
deviatoric  stress 
surface  elevation 
basal  depth 
stress  ratio 


m/s 
kg/m3 
m/s2 
kg/m3 
Pa  s'l/n 


Symbol 


LcA 

L. 

H«A0 

<VUo 

Wk 

*1 


pgH0(i-/>//>j 


W-0)1/nBo 

*i9Ho 

VLo 

VL. 

uyL0 

Bo/miyy1-1'"} 

pgH0 

H 


non-dimensional  pgH0/{B0(U0/L0) 1/n} 


Value 


1 03 

1 .58xl0's 
103 

1 .58xl0's 
104 

6. 33xl08 
104 
103 
10- 1 
6.33x10s 
7.92xl0'9 
917 


1 . 5xl08 


9.0xl06 
1.58x10" 9 
1.58x10- 9 
1.58x10- 8 
5.5xl013 
9.0xl06 
103 
103 

51.4 


Figure  1.  Cross  section  of  ideal  ice  shelf  considered  in  this  study.  The 
trajectory  of  the  ice  column  labeled  £o  follows  one  of  the  two  types  of 
characteristics  associated  with  the  governing  equations.  If  ice-stream 
discharge  was  impulsively  changed  between  two  otherwise  steady  conditions  at 
the  time  this  ice  column  was  at  the  grounding  line,  the  subsequent  locations 
of  this  ice  column  will  separate  the  region  where  steady-state  conditions  have 
been  renewed  (upstream)  from  the  region  where  adjustment  is  still  underway 
(downstream) . 

Figure  2.  If  ice-stream  discharge  (grounding  line  thickness  and  velocity)  is 
changed  impulsively  at  t=0,  unsteady  conditions  will  prevail  only  within  the 
region  of  the  x,t  plane  bounded  by  the  (--constant  and  r*constant 
characteristics. 

Figure  3.  Ice  thickness  (top)  and  flux  (bottom)  anomalies  (non-dimensional 
units,  Table  I)  resulting  from  a  limited  episode  of  modified  ice-stream 
discharge  (scenario  I).  The  ranges  of  x  and  t  displayed  above  are  both  0.0-*1.0 
(non-dimensional  units).  Comparison  with  the  characteristics  shown  in  Figure  4 
suggests  that  thickness  anomalies  primarily  follow  ^-constant  characteristics; 
whereas  flux  anomalies  primarily  follow  r-constant  characteristics. 

Figure  4.  Ice-shelf  response  can  be  organized  in  terms  of  steady  and  unsteady 
behavior  using  the  £-  and  r-characteristics  to  partition  the  x,t  plane.  In 
scenario  I,  for  example,  the  x,t  plane  is  divided  into  five  regions  by  the 
four  characteristics  that  emanate  from  x-0  at  the  start  and  end  of  the  ice- 


stream  discharge  episode.  Regions  I,  III  and  V  are  in  steady  state  with 
current  ice-stream  discharge  (for  region  III,  this  discharge  is  greater  than 
for  regions  I  and  V);  regions  II  and  IV  display  unsteady  conditions. 

Figure  4.  Ice  thickness  anomalies  driven  by  periodic  fluctuations  in  ice- 
stream  thickness  (top,  scenario  II)  and  velocity  (bottom,  scenario  III). 
Contour  intervals  are  0.01  non-dimensional  units,  negative  values  are 
contoured  with  dashed  pattern.  The  ranges  of  x  and  t  displayed  are  both 
0. 0-1.0.  Note  that  maximum  h  for  scenario  III  occurs  downstream  of  the 
grounding  line. 
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